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The two-dimensional Heisenberg model at zero temperature with a quan- 
tum S=l/2 impurity spin, coupled to one site, is studied. The ground state 
properties of the model are calculated for large coupling, in order to study 
the impurity - host spin local singlet formation and related local suppression 
of the magnetization in the host lattice. Analytic results are obtained by de- 
veloping perturbation theory around the exactly solvable two-body (impurity 
plus neighbor) problem and compared with results obtained by exact diago- 
nalization of small clusters. We find that perfect screening of the impurity is 
only achieved for infinitely large interaction, while at intermediate coupling 

the local magnetization is suppressed, but non-zero. 
PACS numbers: 75.10.Jm, 75.30.Hx, 75.50.Ee 



I. INTRODUCTION 



The problem of magnetic impurities interacting with a system of strongly correlated 
electrons has attracted a lot of interest recently, mainly due to the experimental discoveries 
of the high- Tc superconductors and new heavy fermion compounds. In the field of the high-Tc 
materials, the parent compounds are known to be two-dimensional (2D) antiferromagnetic 
(AFM) Mott-Hubbard insulators which are driven to a superconducting state by doping 
(e.g. with holes)0'i. Even though the holes can hop, thus destroying the AFM long-range 
order (LRO) and causing the development of superconducting pairing, the extreme limit of 
static holes is also believed to have physical relevance. More generally, the effect of local 
perturbations on the AFM order is an interesting problem by itself. Several problems in 
2D have been studied - static vacancy (missing site)il, an impurity spin with an on-site@ 
and sublattice symmetricili coupling, as well as an isolated ferromagnetic bondlBi. Impurity 
spins, coupled to ID antiferromagnets have also been considered, by applying bosonization 
techniques!!. 

From the perspective of the heavy fermion physics, in light of some recent experimental 
observations, it is important to take into account strong correlations between the conduction 
electrons, interacting with an on-site magnetic impurity. As a starting point, a model at 
half filling, with a large Hubbard on-site repulsion, leading to the freezing of the conduction 
degrees of freedom, has been proposedi. The problem reduces to an impurity spin, coupled 
to a Heisenberg antiferromagnet. In 2D, for T = the problem has been studied by using 
the linear spin-wave approximation (LSWA)i. For T 7^ 0, when LRO is absent in 2D, the 
Schwinger boson mean-field theory was appliedS. Also, the one-dimensional version of the 
problem has been studied numerically as a toy modellllE. 

In the present work we consider a magnetic impurity, coupled via an on-site Kondo 
coupling to a 2D quantum Heisenberg antiferromagnet at T = 0. Our main goal is the 
investigation of the interplay between LRO and local screening of the impurity spin (Kondo 
effect). Since LRO is present in 2D, with a staggered moment reduced to about 61 % from 
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its classical value, the Kondo coupling, necessary to induce considerable impurity screening 
is rather large. Previous works have used the LSWA with subsequent perturbation theory in 
the Kondo interactioni'0 - an approach not suitable for the study of strong interactions. To 
get insight into the strong-coupling regime, we observe that due to the local character of the 
perturbation it is sufficient first to diagonalize exactly the two-body system, consisting of 
the impurity and directly coupled substrate spin, and then treat the remaining interaction 
with the AFM environment in perturbation theory. Technically we work with a two-particle 
Green's function which we use to construct the ground state of the system. The rest of the 
paper is organized as follows. Our approach is described in Section II. We also compare 
our analytical results with results obtained from exact diagonalization studies. The numer- 
ical procedure is outlined in Section III. Section IV contains summary of our results and 
discussion. 



II. IMPURITY SPIN IN THE HEISENBERG ANTIFERROMAGNET. 

Consider the two-dimensional spin-1/2 Heisenberg model at zero temperature with a 
single, spin-1/2 magnetic impurity: 

H = J J2 ^i-^j + KSo.a (1) 

All the couplings are antiferromagnetic J > 0, A' > 0. We set J = 1 from now on. In the 
first term in (1) the summation is over nearest neighbors on a square lattice. The second 
term represents an S=l/2 spin a coupled to the site i = of the host lattice via the Kondo 
coupling K. It is knownEI that for A' = there is LRO in the ground state of H and the 
staggered magnetic moment m =< >= 0.303. The Kondo interaction favors singlet 
formation and thus suppresses the magnetization locally, leading also to screening of the 
impurity spin. 

In order to treat the strong coupling case A' > 1 we start by rewriting the Hamiltonian 

as 
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(2) 



Hsw = E ^k(ak"k + /^k^k)- (3) 

k 

The spins Si,i = 1 — 4 are the nearest neighbors of 5*0 (see Fig.l.) and are assumed to 
belong to sublattice B (spin down). Here /?k and a^. are the usual spin- wave operators with 
dispersionS: 
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£k = 2y 1 - ll, 7k = -(cos(A;:r) + cos{ky)). (4) 

In (2) we have explicitly separated the host spin-impurity and host spin-AFM background 
interaction. We have assumed that the exclusion of one spin (the one at i = 0) does not 
influence the spin-wave Hamiltonian. This is valid in the one-loop approximation (see below) 
which we use in the present work. However, if one wants to go beyond the lowest order (i.e. 
to higher loops), the influence of the i = perturbation on the spin- wave spectrum has to 
be taken into account. 

Since the impurity spin disturbs the AFM background only locally, we can make the 
mean-field substitution SqS^ Sq < S- >= S^m. The terms in the first curly brackets 
in (2) then can be diagonalized exactly, while the ones in the second curly brackets will be 
treated as a perturbation. For the spins Sq and a we introduce the fermionic representation: 

^o^ = ^!^i, 5o = ^(^{^T-M^i), (5) 

where \E'| and \1'| create (acting on their common vacuum) a Sq component 1/2 and —1/2 
respectively. In order for the two fermions to represent a spin-1/2 operator, they have to 
satisfy the constraint: \E'|\E'| + ^E'j^E'j, = 1, i.e. the physical states have only one fermion. 
The operators representing a via a formula identical to (5) are denoted by $j. Using this 
representation and the usual spin wave expansion for the spins Si,i = 1 — 4, the Hamiltonian 
becomes: 

H = Hsw + Hq + Hint-, (6) 



Ho = 2m(vl>tv[/^ _ ^^^J^) + :^(x[/{v^^ - v[/|v[/^)($|$^ - $1$^) + ^{^^H^^]^^ + h.c), (7) 



H^nt = \ T7^T^i E 7k(Mk/5k + t^kalk) + h-c, (8) 



where N is the total number of lattice sites, and -Uk = ~ ~sign(7k)y — | + ^ 

are the parameters of the Bogoliubov transformation. The k sums in (8) as well as all future 
formulas are over half of the first Brillouin zone, < k^, ky < vr. 

We have previously apphed an approach, similar to the one that leads to the effective 
Hamiltonian (6), to study locally frustrating defects in quantum antiferromagnet^l. The 
following treatment is also closely related to the one used by us in Ref.[6]. 

In order to diagonalize Hq it is convenient to define the two-particle matrix Green's 
function G,, 
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(9) 



As (9) suggests, the diagonal elements (11 and 22) correspond to the two-particle states 
|1 >= I ti> and |2 >= | |t>) respectively, where the first arrow represents the spin Sq and 
the second one - the spin a, or, equivalently | ti>= ^1*^110 >> I iT>= ^I'^jl'-' ^' where |0 > 
is the fermionic vacuum. The off-diagonal components represent transitions between these 
two states. For future purposes we define also the states: |3 >= | ]]> and |4 >= | [[>. 
The unperturbed Green's functions (corresponding to Hq) are: 

Next, we evaluate the self-energy corrections to (10) to lowest, one-loop order in perturbation 
theory with respect to Hint- Evaluating the diagram of Fig. 2a, we obtain: 

Here D^^(k,uj) = u — e^^ is the spin- wave Green's function, and G^^{uj) = u — K/A — 2m 
corresponds to the two-particle state |4 >. Analogous calculation gives (see Fig. 2b): 



and Si2 = S21 = 0. Higher loop corrections to the self-energy can also be taken into account. 
Let us mention that vertex corrections do not exist due to spin conservation (reflected by the 
structure of the interaction (8)). The only remaining diagrams are the "rainbow" one M. We 
find that their contribution is small compared to the dominant one of Eq.(ll) and Eq.(12). 
Thus we restrict ourselves to the one-loop order. 

Since the z component of the total spin is conserved, the wave function of the spins 5*0 
and a is spanned by the states |1 > and |2 >, corresponding to Sq + = 0. The equation 
for the effective energy level e* is: 

£* + K/4 + 2m- Sii(£*) K/2 

= 0. (13) 

K/2 e* + K/A -2m -^22 (e*) 

The correct normalized eigenstate is of the form |12 >= > +/i2|2 >, where {^1,^2) is 
an eigenvector of the matrix in (13), and + ^\ = 1. The ground state wave function of 
the two-particle - spin-wave system can be written as 

\G >= v^|12 > |sw > +^5k| ii> pl\sw > +Y.A^\ TT> 4lsw > . (14) 

k k 

We have defined: 

B = /^i7kMk , ^ /X27k^k 

^ ^ Ne* -K/A-2m-ek ^ \ N e* - K/A + 2m - 

and |sw > is the spin-wave vacuum, i.e. ak|sw >= /9k|sw >= 0. The first term in (14) is 

the coherent part of the wave function while the rest is the contribution of the intermediate 

states. The normalization factor Z is defined as: 

The average of any spin operator in the state \G > can be computed by using the explicit 
form (14). The results for the magnetization of the impurity spin and its neighbor as well 
as their spin-spin correlation function are: 
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M(0) >= i^^Z + ^ E (^k - Bl) , (17) 



M(a) >= ^^^Z + \Y.[aI- Bt) , (18) 

k 

1 — 2Z 

C((T, 0) a.So >= ^ + /ii/i2^. (19) 

In order to find the spin-spin correlation between the host spin 5*0 and its nearest neigh- 
bor 5*1 of the host lattice (see Fig.l), we need to find how 5*1 acts on the spin- wave states 
|sw >, aj^|sw >, /5j!;.|sw >, which appear in \G >. This is accomplished by using the Holstein- 
Primakoff representation = — 1/2 + b\bi, Si = b\, and performing the Bogoliubov trans- 
formation 6k = Mk/3k + "^kC^lk- '^^^ final result is 

C(l, 0) Si.So >=< "^o > {-^ + f E + ^ [/^?v^Sn + /i^v^S22] • (20) 

In Eq.(20) the self-energies are evaluated at e = e*. The expression in the curly brackets is 
the LSWA magnetization m = 0.303. In order to compute (17-20) all the lattice sums as well 
as the solution of Eq.(13) have to be calculated numerically. The results are summarized in 
Figures 3 and 4. 

Notice that in the approximation we have adopted for the Hamiltonian Eq.(6), the mag- 
netization of the spins 1,2,3,4 (see Fig.l) is m = 0.303 (since they are part of the spin-wave 
background), and does not depend on K. In order to calculate, e.g. < 5*^ > more accurately, 
one needs to diagonalize exactly the system of three spins SQ,Si,a and treat all the rest 
of the spins in the linear spin-wave approximation. Conceptually the calculation is very 
similar to the one presented above. However, since the technical details are more involved 
and not particularly illuminating, we do not describe them here0. The results for < Sf > 
are given in Table.I. The diagonalization of the three-body system is expected to produce 
also more accurate results for C(l, 0). We have found, however, that the difference between 
the three-body calculation and Eq.(20) is numerically very small. 
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In our technique the magnetization of the unperturbed Heisenberg antiferromagnet {K = 
0) differs from the LSWA value m = 0.303 and depends on the size of the "cluster" - the 
number of spins that are exactly diagonalized. We find the magnetization to be 0.272 if 
the cluster consists of one spin (M(0) from Fig. 3.) and 0.318 for a cluster of two spins 
(< Sf > (GF) from Table. I). We generally expect that the accuracy of our method will 
increase as the size of the cluster increases and higher orders of perturbation theory in the 
cluster-AFM background interaction are taken into account. 

The results for the magnetizations as functions of K (Fig. 3) are consistent with the 
physically expected behavior. As K increases, the impurity spin becomes gradually screened, 
while the magnetization of the nearest substrate spin decreases. Perfect singlet formation 
between 5*0 and a is achieved only in the limit of infinitely large K. The behavior of the 
correlation functions (Fig.4) supports the above picture. The maximum value of 0.75 for 
the correlation function C{a,0) is gradually approached with the increase of K. Notice 
that < > (Table. I) at any 7^ is larger than the value at = 0. This means that 
the quantum fluctuations have effectively decreased at that site. Such behavior was also 
observed in the 2D Heisenberg model with a missing sitei. 

For small value of the Kondo coupling K < 1 our results for the magnetization are 
consistent with previous calculations of Igarashi, Murayama, and FuldJ in this regime, who 
used LSWA in combination with perturbation theory in K. We emphasize, however, that 
our calculation was specifically designed to treat the strong-coupling limit K > 1. Similar 
results to ours have been reported in the one- dimensional version of the model by Igarashi 
et alM. They studied numerically the behavior of the quantities, defined in Eqs.( 17-20). 
From comparison of our curves with theirs, we conclude that the tendency for local singlet 
formation is more pronounced for a Heisenberg chain, which is expected, since LRO is not 
present in ID. 
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III. EXACT DIAGONALIZATION STUDIES. 



We have also studied the system by exact diagonahzation of small clusters. Traditionally 
periodic boundary conditions (PBC) are usedEl. However we have found in the study of 
related problems, involving locally frustrating spin defect^l that the use of PBC may lead 
to spurious discontinuities in the observables. We have attributed this to finite size effect^ . 
Instead, in our previous worki we applied a staggered magnetic field in the z direction to 
spins on the boundary of the cluster (see Fig.l). The same procedure is used here. We 
have chosen a cluster of = 18 + 1 spins and the boundary field to give < >= 0.3 on 
the boundary spins. The advantage of such an approach is that the sublattice symmetry 
is broken which allows us to compute single spin averages and to distinguish between lon- 
gitudinal and transverse correlations. Our results are summarized in Table. I and Figures 
3 and 4, for comparison with the analytic approach. Due to the asymmetry of the cluster 
the correlation function ^(1, 0) is computed as [C{1, 0) + C{2, 0) + C{3, 0) + C{A, 0)]/4 (see 
Fig.l). Analogous symmetrization was used for < Si >. We find a very good qualitative 
agreement between the analytical and exact diagonahzation results. 

Our numerical procedure also breaks the translational invariance of the cluster and will, 
in general, lead to increased finite size effects and slow convergence to the bulk limit. To 
achieve complete self-consistency, the extrapolation to the thermodynamic limit has to be 
performed, by adjusting simultaneously the boundary field to be equal to the magnetization 
of the spins inside the cluster. This certainly would increase even further the numerical 
agreement between the analytical and numerical results. 

IV. SUMMARY AND DISCUSSION. 

To summarize, we have studied the competition between LRO and local singlet formation 
in the 2D Heisenberg model at T = 0, coupled to a magnetic impurity via an on-site Kondo 
term. We were particularly interested in developing a formalism to treat large Kondo cou- 
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plings, since it is physically clear, that due the presence of LRO in 2D, only large couplings 
can lead to substantial impurity screening and local suppression of the magnetization. We 
find that the local singlet is formed asymptotically. At intermediate couplings the magneti- 
zation sustains a non-zero value, which gradually decreases as the coupling increases. This 
picture is supported by the two methods we have used: 1.) An analytic approach, which 
treats exactly the system, consisting of the impurity and the directly coupled to it host spin, 
and takes into account the interaction with the AFM background perturbatively, and 2.) 
Exact diagonalization of small clusters. 

A related problem is the behavior of a Kondo moment in a system without LRO in the 
ground state (non-zero gap in the magnon energy spectrum). Destruction of LRO can be 
achieved e.g. by doping, or by introduction of additional spin interactions. The lack of LRO 
would lead to different behavior, compared to the one found in this paper, since the impurity 
could become completely screened at intermediate couplings. An interesting question is the 
existence of impurity induced bound states in the gap. Such states have been found in 
the S = 1 Heisenberg chainlll, and in s-wave (as well as d-wave under certain conditions) 
superconductors with magnetic impurities^. We plan to address these issues in future work. 

ACKNOWLEDGMENTS 

We would like to acknowledge the financial support of the Australian Research Council. 



1 See, e.g., E. Manousakis, Rev. Mod. Phys. 63, 1 (1991). 

2 For a review, see E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

^N. Bulut, D. Hone, D. J. Scalapino, and E. Y. Loh, Phys. Rev. Lett. 62, 2192 (1989). 
^N. Nagaosa, Y. Hatsugai, and M. Imada, J. Pys. Soc. Jpn. 58, 978 (1989). 



10 



^ J. Oitmaa, D. Betts, and M. Aydin, Phys. Rev. B 54, 15171 (1995). 
^V.N. Kotov, J. Oitmaa, and O. Sushkov, preprint 3ond-mat/9712281 



^K. Lee and P. Schlottmann, Phys. Rev. B 42, 4226 (1990). 

^ S. Eggert and I. Affleck, Phys. Rev. B 46, 10866 (1992); D. G. Clarke, T. Giamarchi, and B. I. 
Shraiman, Phys. Rev. B 48, 7070 (1993). 

^ J. Igarashi, K. Murayama, and P. Fulde, Phys. Rev. B 52, 15966 (1995). 

^'^ K. Murayama and J. Igarashi, J. Phys. Soc. Jpn. 66, 1157 (1997). 

J. Igarashi, T. Tonegawa, M. Kaburagi, and P. Fulde, Phys. Rev. B 51, 5814 (1995); W. Zhang, 
J. Igarashi, and P. Fulde, Phys. Rev. B 54, 15171 (1996). 

12 W. Zhang, J. Igarashi, and P. Fulde, J. Pys. Soc. Jpn. 66, 1912 (1997). 

1^ The absence of vertex corrections in this problem is analogous to the same effect in the t-J 
modeli, and follows from the form of the fermion - spin-wave vertex. 

1^ For details of a similar calculation, involving three spins, see Ref.[6], where the problem of an 
external spin, coupled symmetrically to the two sublattices, is studied. 

15 J. Oitmaa and D. Betts, Can. J. Phys. B 56, 897 (1978). 

16 E. Sorensen and I. Affleck, Phys. Rev. B 51, 16115 (1995). 

17 Pqj, recent work and references, see M.I. Salkola, A.V. Balatsky, and J.R. Schrieffer, Phys. Rev. 
B 55, 12648 (1997). 

TABLE I. Magnetization on site 1 calculated by the Greeen's function (GF) technique and by 
exact diagnalization (ED) of N=18+l cluster. 
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1 < Sf 


> l(GF) 


0.318 
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0.329 


0.328 


1 < -^i 


>|(ED) 


0.3717 


0.3804 


0.3822 


0.3820 


0.3814 


0.3809 
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FIGURE CAPTIONS 

FIG. 1. Schematic representation of the system. The black dots are the host spins, with 
interaction J=l between them. A cell, containing N= 18 + 1 sites, used for exact diagonalization, 
is shown. 

FIG. 2. One-loop diagrams, contributing to the two-particle Green's function, Eq.(9). Solid 
and dashed lines represent fermion and spin-wave Green's functions, respectively. 

FIG. 3. The impurity spin average M{a) and the magnetization of the neighboring host spin 
M(0). The solid and dashed curves are calculated from Eq.(17) and Eq.(18), while the circles and 
squares are the corresponding exact diagonalization results. 

FIG. 4. Spin-spin correlation functions as defined by Eq.(19) and Eq.(20). The solid and dashed 
line are obtained by numerical evaluation of Eqs. (19-20). The open symbols are the corresponding 
exact diagonalization results. 
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